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Motivation 


Introduction and Motivation 
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Motivation 


Path integral with complex weights appear in many important 
physics: 


ə Finite-density lattice QCD, 
spin-imbalanced nonrelativistic fermions 


ə Gauge theories with topological 0 terms 


ə Real-time quantum mechanics 


Oscillatory nature hides many important properties of partition 
functions. 
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Example: Airy integral 
Let's consider a one-dimensional oscillatory integration: 


a= | i T as 
p 3 - 


RHS is well defined only if Ima = 0, though Ai(z) is holomorphic. 


Horo nap 


exp( osx ested) 


AAACN 


Figure : Real parts of integrands fora = 1 (x10) & a = 0.5i 
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Contents 


ə How can we circumvent such oscillatory integrations? 
=> Lefschetz-thimble integrations 
[Witten, arXiv:1001.2933, 1009.6032] 
[YT, Koike, Ann. Physics 351 (2014) 250] 


oe Applications of this new technique for path integrals 


> Study on phase transitions of matrix models. 

> Lefschetz-thimble method elucidates Lee-Yang zeros. 
[Kanazawa, YT, JHEP 1503 (2015) 044] 

> General theorem ensuring that Z € R. 
Sign problem in MFA can be solved. 

> Application of the theorem to the SU(3) matrix model 


[YT, Nishimura, Kashiwa, in preparation] 
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Lefschetz-thimble integrations 


Introduction to Lefschetz-thimble integrations 
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Again Airy integral 


Airy integral: 


3 
Ai(a) = f T epi (5 + ar) 


The integrand is holomorphic w.r.t x 
=> The contour can be deformed continuously without changing the 
value of the integration! 
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Lefschetz-thimble integrations 


Steepest descent contours 
What is the most appropriate contour for our purpose? 


Re[S(x,a)] should be made as small as possible. 
<> The contour should be perpendicular to Re|S(x,a)] = cosnt. 


Steepest descent ones J must satisfy Im[S(x, a)] = const. 
because of the holomorphy. 


m 


[i 
-t 1 P | | 
Figure : Contour plots for Re[S(z, a)] with a= = exp: +i0. 1. 
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Lefschetz-thimble path integral 


Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (o = 1,2) for the Airy 


integral: 
d 3 
Ai(a) = > No [ Pi (3 taz) 


Ng: intersection number of the steepest ascent contour Ke and R. 
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Figure : Lefschetz thimbles J and duals K (a = exp(0.1i), exp(7i)) 
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Tips: Airy integral & Airy equation 


Let us consider the “equation of motion”: 


dz d ei(2?/3+a2) 


= 0. 
Qn a 


This is nothing but the Airy equation: 


Two possible integration contours 
= Two linearly independent solutions of eom. 
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Lefschetz-thimble integrations 


Generalization to multiple integrals 


Model integral: 
z= | dzı -- - d£n exp S(x). 


What properties are required for Lefschetz thimbles J? 


@ J should be an n-dimensional object in C”. 
© Im[S] should be constant on J. 
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Lefschetz-thimble integrations 


Short note on technical aspects 


Complexified variables (a = 1,..., n): Za = Za + ipa- 


Regard x, as coordinates and p, as momenta, 
so that Poisson bracket is given by 


N a 


a=1,2 


O1 Opa Oxa Opal 
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Lefschetz-thimble integrations 


Short note on technical aspects 


Hamilton equation with the Hamiltonian H = Im[S(z.)]: 


df (x, p) dza as 
Ta SI (= dt --(#)) 


This is Morse’s flow equation (= gradient flow): 


(=) 


a= 


<0 
dt 


= We can find n good directions for J around saddle points! 
[Witten, 2010] 
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Multiple integral: Lefschetz-thimble method 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles Js: 


d"r & = ne | d”z eS), 
h. ee), 


Js are called Lefschetz thimbles, and Im[S] is constant on it. 


Ne. intersection numbers of duals X, and R”. 
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Gross-Neveu-like model 


Phase transition associated with symmetry 
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Gross-Neveu-like model 


0-dim. Gross-Neveu-like model 


The partition function of our model study is the following: 


N 


Zu(G,m) = f addy exp IT tmt ( Sia)’ | Ä 


a=1 
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Gross-Neveu-like model 


Hubbard-Stratonovich transformation 


N 
Zy(G,m) = 4/ “a a 


G 


For simplicity, we put M = O in the following. 


Bosonization: 


with 


S(o) log[p? +(e + m)?]. 


S has three saddle points: 


_ OS(z) _ 2z 2z 


0 = => 2g=0, +VG-P. 


dz G p+ 
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Behaviors of the flow 
Figures for G = 0.7e-° 4, 1.1e7”! at p? = 1 
[Kanazawa, YT, arXiv:1412.2802]: 


ə z = 0 is the unique critical point contributing to Z if G < p? 
ə All three critical points contribute to Z if G > p°. 
Yuya Tanizaki (University of Tokyo, RIKEN) 


Lefschetz-thimble path integral 


DAG 


Stokes phenomenon 


The difference of the way of contribution can be described by 
Stokes phenomenon. [Witten, arXiv:1001.2933, 1009.6032] 


Figures of G-plane for ImS(0) = ImS (z4) [Kanazawa, YT, arXiv:1412.2802]: 


-1.5 


=1,0:=0:5: 0.0 05 10 15 20 


RG 
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Gross-Neveu-like model 


Dominance of contribution 


The Stokes phenomenon tells us the number of Lefschetz thimbles 
contributing to Zy. 
However, it does not tell which thimbles give main contribution. 


Zn ~ #exp(-NS(0)) + #exp(—NS(z4)) 


In order to obtain (a) Æ 0 in the large-N limit, z should dominate 
z=0. 


> ReS(z; 


) < ReS(0) 
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Connection with Lee-Yang zero 


Blue line: ImS(z+) = ImS(0). 
Green line: ReS(z;) = ReS(0). 
Red points: Lee-Yang zeros at N = 40. [Kanazawa, YT, arXiv:1412.2802] 
Apr. 3, 2015 @ YITP, Kyoto 21/39 


Gross-Neveu-like model 


Conclusions for studies with GN-like models 


@ Decomposition of the integration path in terms of Lefschetz 
thimbles is useful to visualize different phases. 


@ The possible link between Lefschetz-thimble decomposition and 
Lee-Yang zeros is indicated. 
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Sign problem in MFA 


Sign problem in MFA and its solution: Formal discussion 
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Sign problem in MFA Formal discussion 


Airy integral with real a 


If a € R, the following partition function is real: 


1 x? 
Zairyla) = z [ ropi (= + ar) A 


However, the integrand is complex-valued. 


Why does this happen? 
What is the mean-field approximation of this partition function? 
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Sign problem in MFA Formal discussion 


Airy integral: Antilinear symmetry 


The system have an antilinear symmetry: 


By summing up in pairs +2, the result becomes manifestly real! 
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Airy integral: Lefschetz thimble 


Extend the linear map x ++ —x to the antilinear map z +> — 2. 


The Lefschetz thimbles form invariants under the antiliner map: 


M 
Figure : Lefschetz thimbles 7 and duals K (a = exp(0.1i), exp(7i)) 


ak! 
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Sign problem in MFA Formal discussion 


Generalization: Set up 


Consider the oscillatory multiple integration, 


z= | d”rexp-S(z). 


To ensure Z € R, suppose the existence of a linear map L, satisfying 


ə S(x) = S(L- x). 
oi? =1, 


Let’s construct a systematic computational scheme of Z with 
ZER. 
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Sign problem in MFA Formal discussion 


Generalization: Flow eq. and Complex conjugation 


Morse’s downward flow: 


Complex conjugation of the flow: 


dt Oz; 


and thus z’ := L -Z satisfies Morse’s flow equation . 
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Sign problem in MFA Formal discussion 


Generalization: Immediate conclusion 


Denote the antilinear map as K(z) := L-Z. 
@ When z° is a saddle point, so is K(z7). 
@ Let J, is the Lefschetz thimble around 2°. 


=> JE :={K(z) | z € Je} coincides with the Lefschetz thimble 
around K(z”) (up to orientation). 


© Intersection numbers are the same: (Ko, R”} = (KE, R”) under 
appropriate choice of the orientation. 


@ The contributing thimbles always form an invariant pair 


ge. 
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Sign problem in MFA Formal discussion 


General Theorem 


Let us decompose the set of saddle points into three parts, 
= UN}; UX, where 


Xo = {o| = L-27}, 
{o | ImS(27) > 0}, 
x. = {o |ImS(2z7)< 0}. 


M 
+ 
| 


The antilinear map gives X, ~ %_. 
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General Theorem 


The partition function: 


Z = ne | d”z exp —S(z) 
oEX0 To 
+ n f d"z exp —S(z) 
fen, TIA 


Each term on the r.h.s. is real. 


(YT, Nishimura, Kashiwa, in preparation) 
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Sign problem in MFA 


Sign problem in MFA and its solution: 
Application to QCD-like theories 
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Aue 2 
The QCD partition function at finite density 


The QCD partition function: 
Zacn = [ra detM (tax; A) exp —Sym [A], 


w./ the Yang-Mills action Sym = 3tr i da? f d°x|F,,|? (> 0), and 


detM(u, A) = det |y” (3, + igAL) + Vlak + Max] - 


is the quark determinant. 
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SW HELEN TE Application to QCD 


Charge conjugation 


If Hax # 0, the quark determinant takes complex values 
=> Sign problem of QCD. 


However, Zoco € R is ensured thanks to the charge conjugation 
Ar —At: 


det M(jigk, A) = det M(—jigx, A") 
det M (fax, — A). 


(First equality: y5-transformation, 
Second equality: charge conjugation) 
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SruweelJlumumWii\g Application to QCD 


Polyakov-loop effective model 


The Polyakov line L: 


l 1 i it i l i 
L= z dag [e@1+92), eil 01t0a) p20], 
Let us consider the SU (3) matrix model: 
m [wono 02) exp [Ver (81, 82)] , 


where H = sin? 0, sin? (04 + 305) /2 sin? (01 — 302) /2. 
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Smee Vigame Application to QCD 


Charge conjugation in the Polyakov loop model 


Charge conjugation acts as L + Li: 
Vor (21, 22) Zu Ver(ZT, —%2). 


Simple model for dense quarks (£ := trL): 


37-1 
Var = HET (erta(01,02) + etal, 0) 
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SW HELEN TE Application to QCD 


Behaviors of the flow 


Ot peo kee ee oe He ee 


O 2 


7 
Re(zı) 
Figure : Flow at h = 0.1 and u = 2 in the Re(z1)-Im(z>) plane. 
The black blob: a saddle point. 
The red solid line: Lefschetz thimble 7. 
The green dashed line: its dual K. (YT, Nishimura, Kashiwa, in preparation) 
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SruwelelJlumlm Wi Application to QCD 


Saddle point approximation at finite density 


The saddle point approximation can now be performed. 


Polyakov-loop phases (z1, z2) takes complex values, so that 


(4), (0) € 


Since Im(z2) < 0 and 
L~ 


we can confirm that 


(YT, Nishimura, Kashiwa, in preparation) 


R. 
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(2e? cos 6; + eo), 
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Oiane gee 
Summary for the sign problem in MFA 


ə Lefschetz-thimble integral is a useful tool to treat multiple 
integrals. 


ə Saddle point approximation can be applied without violating 
ZER. 


ə Sign problem of effective models of QCD is (partly) explored. 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Apr. 3, 2015 @ YITP, Kyoto 39 / 39 


